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Abstract 
It  is  shown  that  macroscopic  correlations  in  a  fluid  are 
conserved  for  macroscopically  long  times.   The  equations  of 
conservation  can  be  written  in  a  form  independent  of  the 
density  of  the  fluid  and  are  therefore  valid  for  a  liquid  as 
well  as  a  gas.   For  the  purpose  of  describing  correlations  in 
a  rarefied  gas,  a  generalization  of  the  Boltzmann  equation  is 
suggested  which  not  only  gives  the  correct  macroscopic  corre- 
lation equations  but  satisfies  certain  compatibility  criteria. 
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I.    Introduction 

The  aim  of  kinetic  theory  and  statistical  mechanics  is 
to  explain  the  macroscopic  properties  of  matter  on  the  basis 
of  the  properties  of  the  (microscopically  small)  molecules  of 
which  it  consists.   Accordingly,  in  kinetic  theory,  a  fluid 
is  regarded  as  a  dynamical  system  consisting  of  a  very  large 
number  of  particles  —  we  shall  use  the  terms  "molecule"  and 
"particle"  synonymously.   Two  molecules  of  a  fluid  are  said 
to  be  correlated  if  the  position  and  velocity  of  one  of  them 
are  not  independent  of  those  of  the  other. 

Correlations  are  of  two  kinds  ~  short-range  and  longe- 
range.   When  the  distance  of  separation  between  any  two  mole- 
cules is  of  the  order  of  the  range  of  the  intermolecular  force, 
the  correlation  is  said  to  be  short-range;  when  it  is  consider- 
ably larger,  the  correlation  is  long-range.   Short-range  corre- 
lations are  studied  in  the  theory  of  dense  gases;  long-range 
correlations  usually  appear  in  a  context  apparently  unrelated  to 
their  microscopic  nature,  namely,  in  turbulent  fluid  dynamics. 
One  of  our  aims  is  to  clarify  this  connection. 

Contributions  to  the  long-range  correlation  between  any 
two  molecules  come  from  their  initial  correlation,  their  present 
interaction  with  other  molecules,  and  any  collisions  that  may 
have  taken  place  betv;een  them  and  other  molecules  in  the  past. 
The  Boltzmann  equation  (see  Sec.  7  or  Ref.  [1]  for  a  derivation), 
which  describes  a  rarefied  gas,  is  based  on  the  assujnption  that 
long-range  correlations  are  negligible.   (It  is  even  stated  some- 
times  that,  over  m.acroscopic  times,  they  are  destroyed  by 
collisions  [2].)   The  question  naturally  arises  whether  this 
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assumption  is  correct.   It  v-as  conjectured  by  Grad  [3]  that  it 
may  not  always  be,  and  vie   verify  this  conjecture.   More  pre- 
cisely, his  conjecture  is  that  long-range  correlations  are 
conserved  in  the  sense  that  when  they  are  expressed  as  macro- 
scopic correlations,  i.e.,  as  correlations  between  any  two 
conserved  quantities,  they  obey  conservation  equations.   Thus, 
if  the  macroscopic  correlations  are  large  initially,  they  will 

■X- 

continue  to  be  large  for  macroscopically  long  times. 

Long-range  microscopic  correlations  can  be  transformed 
into  correlations  between  macroscopic  quantities  as  follows. 
By  introducing  a  probability  density  on  the  phase  space  of 
the  system  and  by  localizing,  with  respect  to  the  molecules, 
the  macroscopic  quantities  of  interest  to  us  —  such  as  mass, 
momentum  and  energy,  we  can  calculate  their  densities  in 
ordinary  3-space.   Similarly,  given  two  points  in  3-space,  we 
can  find  the  joint  density  at  these  points  of  any  two  (not 
necessarily  distinct)  macroscopic  quantities.   Subtracting 
from  this  joint  density  the  product  of  the  densities  of  the 
two  quantities  at  the  given  points,  we  get  the  correlation 
between  them.   Our  result  is  that,  if  each  of  the  two  m.acro- 
scopic  quantities  is  conserved,  and  if  the  distance  bet\>;een 
the  two  points  chosen  is  larger  than  a  few  times  the  range  of 
the  intermolecular  force,  then  the  correlation  is  conserved. 
In  other  words,  certain  moments  of  long-range  correlations   are 


Correlations,  which  we  sha3.1  study  here,  refer  to  one  time  and 
two  or  more  positions,  v^hereas  fluctuations,  vihich  we  shall  not 
study  here,  refer  to  two  or  more  times  and  any  number  of  positions. 
Fluctuations,  like  correlations,  can  be  studied  either  from  the 
point  of  viev;  of  fluid  dynamics  [4]  or  from  that  of  kinetic  theory 
[5]. 


conserved. 

The  macroscopic  correlation  equations  for  a  rarefied  gas 
have  been  given  by  Zhlgulev  [6]  and  Tsuge  [7,8].   Tsuge  also 
raised  the  question  of  the  validity  of  the  Boltzmann  equation 
in  the  presence  of  turbulence.   It  is  clear  that  the  equation 
is  valid  as  long  as  macroscopic  correlations  are  small,  i.e., 
as  long  as  the  motion  of  the  fluid  is  not  turbulent.   When 
there  is  turbulence,  all  the  derivations  of  the  equation  break 
down,  but  it  is  not  clear  whether  the  equation  Itself  remains 
valid.   At  any  rate,  in  order  to  obtain  the  correct  macroscopic 
correlation  equations  for  a  rarefied  gas,  one  should  replace  the 
Boltzmann  equation  by  a  kinetic   equation  or  sequence  of  equations 
which  Includes  correlations.   This  is  precisely  what  we  do  in 
Section  7.   Our  work  could  thus  be  interpreted  as  a  further 
step  tov;ard  a  kinetic  theory  of  turbulence  along  the  lines  of 
Zhlgulev  and  Tsuge. 

A  few  remarks   are  now  in  order.   Conventional  thermo- 
dynamics is  based  on  the  assumption  that  the  system  has  only 
two  time-independent  integrals,  namely,  mass  and  energy.   If 
the  system  has  other  Integrals,  the  thermodynamics  that  ensues 
differs  from  conventional  thermodynamics  and  has  been  discussed 
by  Grad  [10,11]  and  Lev/is  [12].   Fluid  dynamics,  however,  is 
based  on  the  existence  of  another  Integral,  namely,  momentum. 
Systems  whose  angular  momentum  is  conserved  are  also  sometimes 
considered'-   '   »  Jj   J  ^   pq^.  -^j^^q  sake  of  generality,  we 


The  idea  that  microscopic  correlation  functions  can  be  used  to 
describe  turbulence  occurs  in  Ref.  [9],  but  no  equations  for 
macroscopic  correlations  are  derived. 
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consider  a  system  whose  molecules  have  internal  degrees  of 
freedom  and  v/hich  has  any  or  all  of  the  above  types  of 
integrals.   Also,  our  derivation  is  valid  even  when  the  inter- 
molecular  forces  are  velocity  dependent,  though  we  require 
that  the  divergence  v;ith  respect  to  velocity  of  the  force  vanish. 

The  method  we  use  is  essentially  the  same  as  that  developed 
by  Grad  [11]  in  his  generalization  of  Irving  and  KirlciA/ood ' s  [15] 
derivation  of  the  conservation  equations  of  fluid  dynamics. 
Grad's  method  is  more  general  than  that  of  Irving  and  Kirkwood 
in  two  ways:   (1)  Instead  of  considering  only  mass,  momentum, 
and  energy  as  Irving  and  Kirkwood  did,  he  considered  general 
integrals  v;hich  are  purely  additive  —  as  mass  and  momentum 
are  —  or  which  have  a  purely  additive  part  and  a  part  that 
consists  of  two-body  interactions,  as  energy  does.   (2)  He 
considered  systems  v/hich  have  internal  degrees  of  freedom. 

Because  of  the  generality,   Grad's  method  is  readily 
applicable  to  our  problem.   Essentially,  it  consists  in 
computing  the  density,  in  3-space,  of  an  arbitrary  integral, 
multiplj'-ing  the  Liouville  equation  by  it,  and  integrating  the 
resulting  equation  with  respect  to  everything  except  the  space 
coordinate.   What  one  gets  then  is  an  equation  for  the  rate 
of  change  of  the  density  of  the  integral.   In  order  to  show  that 
it  is  a  conservation  equation,  one  has  to  just  show  that  the 
terms  other  than  the  convectlve  derivative  of  the  density  are 
divergences  of  some  vector  functions.   To  derive  the  conservation 
equations  for  macroscopic  correlations,  we  proceed  in  a  similar 
manner,  lie   first  compute  the  joint  density  in  3-space  of  a 
pair  of  integi-als,  then  by  multiplying  the  Liouville  equatj.on 


by  it  and  integrating  with  respect  to  everything  except  the 
two  space  coordinates  involved,  obtain  an  equation  for  the 
rate  of  change  of  the  density.   It  turns  out  that,  if  the  two 
space  coordinates  are  separated  by  a  sufficiently  large  distance, 
the  rate  equation  becomes  a  conservation  equation.   Using  this 
equation  and  the  equations  for  the  conservation  of  the  density 
of  each  of  the  two  Integrals,  we  get  the  conservation  equation 
for  the  correlation  between  them. 

We  assume  that  the  molecules  of  the  fluid  obey  the  laws  of 
classical  mechanics.   Also,  for  the  sake  of  completeness  of 
presentation,  we  will  reproduce  the  relevant  portions  of 
Grad's  paper.   In  particular,  the  material  of  Sections  4  and 
5  is  taken  from  his  paper. 

I  take  this  opportunity  to  express  my  sincere  gratitude 
to  my  teacher  Prof.  Harold  Grad,  who  not  only  suggested  the 
problem  but  provided  constant  encouragement  and  guidance. 
Also,  I  thank  Miss  Rachelle  Rousso  for  her  painstaking  typing 
of  the  manuscript. 
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2.   Notation  and  Conventions 

We  consider  a  system  consisting  of  N  particles,  each  of 
which  has  s  degrees  of  freedom  —  3  external  and  s-3  internal, 
We  denote  the  coordinates  of  the  i-th  molecule  by 


^.  =  (x. ,  ^. ,  q .  ,p. ) 


where  x.  is  the  position  vector  of  the  center  of  mass  of  the 

^         dx. 

molecule  and  ^.  =  — tt-  is  the  velocity  vector  of  the  center  of 

dq 

mass.   The  q.  denote  internal  coordinates  and  the  p.  =  -— rr- 
^1  ^1    dt 

denote  internal  velocities. 
Let 


P  —  (  r^,  r2)  Po>  •  •  •  J  Pj^j 


denote  a  point  in  the  2Ns-dimensional  phase  space  of  the 
system.   We  introduce  a  probability  density  f(P,  t)  on  this 
space  and  normalize  it: 


/f(P,t) 


dP  =  1  , 


where  the  integration  extends  over  the  entire  phase  space. 
Whenever  we  suppress  the  time  coordinate  -  as  vie   did  above  - 
we  mean  that  all  the  quantities  appearing  are  considered  at 
one  instant  of  time,  viz.,  t.   Nov;,  we  consider  a  one-phase, 
one-component  system,  so  that  all  the  molecules  are  identical. 
So  v/e  assume  that  f  is  symmetric  in  all  the  particles,  and 
define  the  reduced  or  marginal  densities  -  called  the  one-particle 
distribution  function,  the  tv/o-particle  distribution  function, 
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etc.  -  by 


^12(^1'^2^  =  /f  dP2dP2^...dPj^  , 


and  so  on.   Since  we  will  need  only  the  first  four  or  five 
distribution  functions,  we  adopt  the  convention  that  if  r 
is  the  total  number  of  subscripts  of  a  distribution  fimction, 
then  it  is  the  r-particle  distribution  function;  the  subscripts 
themselves  indicate  the  arguments  of  the  function.   For  instance, 
fp2ip-  is  the  three-particle  distribution  function  of  the  particles 
labeled  2,'^/and  5.   Clearly,  the  reduced  densities  are  symmetric 
and  normalized  to  1. 

The  density  f  changes  in  time  but  in  such  a  way  that  it 
satisfies  the  Liouville  equation: 

^+|^(^  •  ^1  +-^  •  Fi  +^  •  Pi  +^  .  G. )  -  0  ,        (2-1) 

where  F.  =  — rr  is  the  acceleration  of  the  center  of  mass,  and 
A^'^  dt 

dp. 

G.  =  — tjt  is  the  internal  acceleration  of  the  i-th  molecule. 
Here,  the  •  refers  to  summation  over  the  three  external,  or  the 
s-3  internal,  degrees  of  freedom. 
We  write 


N  N 

F.  =J      F.  .  ,     G.  =y      ^G.  .  ,  (2-2) 


J^^J       -    J=l 


-o 


where  F. .  is  the  acceleration  of  the  center  of  mass  of  the  i-th 
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molecule  due  to  the  J-th  molecule;  similarly  for  G.  ..   F. .  and 
G. .  are  the  accelerations  of  the  i-th  particle  due  to  external 
forces.   ¥e  assume  that 


F^,  =F,i(P.),    G..  =G^^(P,)  (2-3) 


and 


F.J.  =F^2(Pi'Po),    ^IJ  =^12(Pi'^j)  • 

We  consider  two  cases;  an  integral  e  can  be  of  type  (a) : 

N 

e  =  XZ   ^i  '    e.  =  e-,(P.) 
i=l  ^      ^    ^     ^ 

or  type  (b) : 

N       N 
e  =  Z:  ^i  +IZ^  B.J.  ,    e.  =  e^(P.),  and  e.^  =  e^^iF^,I>.) 

i<J 
We  assume  that  e. .  =  s...   If  V  is  any  phase  function,  i.e., 
xf  f  =   -^(P  ,  Pp,  .  . . ,  P  ) ,  we  define  tv;o  mean  values: 

-         -   dP,  dP^ . . .  dP„ 

?(X,)  =N/^f^-gf- N 

and 

dP.dP^. . .dP 


(.(x^.Xg)  =N(N-l)/;f-~2; 


dXp 


Here,  the  convention  is  that  the  integration  is  performed 
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with  respect  to  all  the  variables  except  those  appearing  in 
the  denominator.   The  expected  value  <^>  of  ^P   is  defined  by 


<^>  =  / 


^fdP  . 
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3.   Densities 

In  order  to  calculate  the  amount  of  s  in  any  region  D, 
we  have  to  localize  it.   This  we  do  as  follows:   e.  is  localized 
at  X . ,  the  center  of  the  molecule,  and  e . .  is  localized  half  at 

J-  X  J 

X.  and  half  at  x  .. 

Let  D, ,  Dp  be  any  regions  in  3-space  and  let  *.. ,  <t>p  denote 

their  characteristic  functions  respectively.   The  number  v   of 

N    ^1 


molecules  whose  centers  lie  in  D,  is  given  by  v   =  >   *  (x.). 

^  ^1        i=l  ^     ^ 

Its  expected  value  is 


<VD^>  =  g^/f*l(x,)dP  =  N/^   fl(Pl)dPi 


X- eD, 


Hence,  the  number  density  v  in  3-space  is  N/f  -^ —  .   The 
number  of  pairs  of  molecules  such  that  one  is  in  D,  and  the 


other  is  in  Dp  is 


JZ  <J'i(x.)*2(x.)  =ZZ  *i(Xi)*2(x.)  +JZ  ^(Xi)<J>2(xJ  . 

The  first  sum  on  the  right  hand  side  takes  into  s,ccount  the 
situation  where  one  particle  is  in  D^  f)  Dp  and  is  counted  as 
a  (spurious)  pair.   By  calculating  the  expected  number  of 
pairs,  we  find  that  the  pair  density  is  given  by 


dP^        f^        dP^  dP^ 


'(^1-^2)K  d5^-'^(^-^)ri2  d5^ 


dXg 


To  calculate  the  density  of  an  integral  e,  suppose  first  that 

it  is  of  type  (a).   Then  the  amount  e_^  of  e  in  D^  is  ^       '^^{x.)e., 

BJid 
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Vl 


Hence,  the  density  of  e  is 


If  e  is  of  type  (b). 


=  n/   e^f^dP^+M|^/   e^gf^^^P^dP^  . 
X,  eD-  x^  eD- 

So,  the  density  of  e  is 

^7  1  1  dx^  ^   2    J^12^12   dx^  ^^  '^^ 

Now,  we  shall  calculate  the  joint  density  of  any  pair  of  integrals . 
There  are  three  cases: 

(1)  Both  integrals  are  of  type  (a). 

(2)  One  integral  is  of  type  (a),  and  the  other,  of  type  (0). 

(3)  Both  integrals  are  of  type  (b). 

Strictly  speaking,  case  (3)  subsumes  cases  (1)  and  (2).   Never- 
theless, for  the  sake  of  convenience  and  clarity,  we  shall  do 
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them  separately. 

Case    (1):  e',    e"   both  of  type    (a). 


^h^=IZ\(\)n  ■ 


=  ^f  4^lVPl  -^N(N-l)  /       4^2'^2dPl^P2    • 

So,  the  Joint  density  of  e'  and  e"  is 

/*      dP^        /*       dP,dP^ 
N5(x^-X2)/e{ 4'fi  ^  +N(N-l)/s' e^f^2  d^  '         (3-3) 

Case  (2) ;      e'  of  type  (a)  and  e"  of  type  (b) . 

°2     j    2   J   J   2  ^   JK  d   J 
<^i^eS^>  =  ^/BlB;:0,(x,)*2(x,)fdP+g:/B.B':.,(x.)0g(x.)fdP 

+  i.^^,/%k^(-i)^2(-j)f^P 

=  N  r      e{e"f-,^dP^  +N(N-l)y    ^{^2^\2'^'^1^^2 

x^eD^ 
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2         J  ^1^12^12^^1^^2  ^^   ^J    ^1^12^12   l^^S 


XgEDg 


^  "'"-^^'^'-^^  /   eis^3f,,3.P,CP,.P. 


X-.  eD^ 


Hence,  the  Joint  density  of  s'  and  e"  is 

dP.         -,       dP^  dP^ 
m(x^'X^)  /e-!e:f.  ttA  +N(N-1)  leJie"f.^  ^"^^^ 


+  M.5=ll  5rx  -X  ^A'e"  f   ^^1^^2+N(N-1)  T,,,..  f   ^^1^^2   (3„4^ 
^^    o^x^  ^2''J^1^12^12   dx^     2    7^1^12^12  dx^dx^  ^^  ^> 

.    N(N-l)(N-2)  r  ,  ..  ^    dP.^_dPgdP^ 

2     ,/^l  23  123   dx^dxp   • 


Case  (3):      e',  e"  both  of  type  (b) . 

Proceeding  as  in  cases  (l)  and  (2),  we  find  that  the  joint 

density  of  e'  and  e"  is 


,(x^-X2)/B{s;;f^  ^+N(N-l)/s'£''f^2    d^ 


2 

dP,  dP^     „/,,  ,v  /%  dP.dP 


N(N-ll,  .r,„  ^^1^.^2      N(N-l)  r.    „  ^^1 

i-         2       o^x^  X2;ye^E-L2^12        dx^     ^      2       J  ^1^12^12   dx^ 


2 


dXg 


dP.dP^dP^     ,,/,,  .^  r^  dP.dPg 


+  MihiJiiiisi r    ,,         ^^i^^2^^3  |N(N-i)g,        s  r     „     ^^^ 

^  2  ./ ^1^23^123        dx^dXg     ^^      ^"^"^1  ^2^./ ^12^1^12      dx^ 


dP.dP^     ^,r.r  nx/.T  ON  /*  dP, dP^dP^ 


2       7^12^2^12   dx^dx^  2  ,/ ^13^2^123     dx^ 

+  — ^Tf       ^°^^i   ^2:/ ^12^12^12        dx^ 


dXg 
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+  _.    Zf"    ^""^^^l  ^2'',/  12  13  123     dx^ 


■Jir^^    ^"  f   ^V^2  ,N(N-l)(N-2)r.   „      ^^1^^^ 
J  ^12^12^12  dx-^dXg        4    7/^13  21^123   dx^c 

^r    ^^/ ^12^23  123   dx^dXg 


,  N(N-lUN-2)  r 


dP^  dP^dP„ 

e"  f    1 2 3 

13  23  123   dx^dXg 

dP.  dP^dP^dP, 


^N(N-1)(N-2)(N..3)/;.  ,n  f   ^  ar^a^gar  ar^ 

^^      4       :/^13  2^  123^    dx^dXg  ''^  ^^ 

The  joint  densities  consist  of  three  parts.   The  terras 
containing  the  5-function  give  rise  to  self- fluctuations  or 
auto-correlations.   Those  in  which  the  molecules  at  x,  and  Xp 
interact  with  each  other  either  directly  or  via  a  third  molecule 
express  short-range  correlations.   The  rest  of  the  terms  refer 
to  long-range  correlations.   Since  we  are  interested  only  in 
long-range  correlations,  we  assume  that  the  molecules  labeled  1 
and  2  are  far  apart.   More  specifically,  if  a   denotes  the  range 
of  the  inter-molecular  force,  we  assume  that  in  case  (1), 
particles  1  and  2  are  separated  by  a  distance  larger  than  a. 
Then  the  density  (3~3)  becomes 


N 


r       dP  dP 
(N-l)/sis^f,2  ^^  (3-6) 


Similarly,  if  we  assume  that  particles  1  and  2  are 
separated  by  a  distance  greater  than  2a  in  case  (2)  and  3^   in 
case  (3),  the  densities  {3-^)    and  (3-5)  become 
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N 


(N-l)/e-' 


p"f   .^^1^2  ,  N(N-lUN-2) 
^2^2  dx^dx"  """"^    2 


:h^ 


'  f. 


dP^dP^dP^ 


23  123   dx^dx^ 


(3-7) 


and 


dP,  dP. 


N 


(''-)Ax4'rx.^-^^i^^^4^/(^ie^,^--  -     '-^^^^ 


'23^2^13^^123   dx^dx~ 


^N(N-l)(N-2)(N-3)/;,  ^„  . 
4 7  13^24^ 


dP^dP^dP^dP^^ 
1234    dx^dXp 


(3-8) 


respectively. 
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4.    Some  Identities  and  a  Formula  of  Irving  and  Kirkwood 

We  shall  now  use  the  fact  that  s  is  an  integral.   First, 
we  introduce  some  notation: 

be.  be. 


and 


be.  .  be. 

ik 


(i,a,k)  =^   .  F.^+^  .  G.. 


1 

there  is  no  summation  on  repeated  indices. 

If  there  are  no  external  forces  acting  on  a  system,  we  say  it 
is  isolated.   Suppose  that  the  system  is  isolated  and  consists 
of  only  one  particle.   Then 


d£,    be  be,  be^  be^ 

dt   dx,    ^1   df,     1  oq  ^1     dp^    1 


.  G,  =  0 


or 


be-.  be.. 

OX-    ^1  bq_j^        ^1  ^    ' 


This  is  an  identity  that  limits  the  class  of  functions  e, 
which  make  e  an  integral.   If  the  system  consists  of  two  par- 
ticles and  is  isolated,  then 


-TT-(  e-,  +£0)  =0      in  case  (a) 


and 
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•^(  e-j_  +  Sg  +  £^2)  =0   Incase  (b)  . 


Using  the  identity  (4-1),  we  get 

(1,2)  +(2,1)  =  0  (4-2) 

in  case  (a),  and 
(1,2)  +{2,l)+(l,2,2)+(2,l,l)+4|i  .  e,+4i|-?2+^-Pl 

as,, 

in  case  (b) . 

Now  suppose  that  the  system  consists  of  three  particles 
and  is  isolated.   Consider  case  (b) .   Then 


j^i  £^+£2+62+  e-L2  ^  ^13  ^  ^23^  ^   ° 


Using  the  identities  (4-1)  and  (4-3),  we  obtain  from  this  the 
identity 

(1,2,3)  +(2,1,3)  +(2,3,1)  +  (3,2,1)  +(3,1,2)  +(1,3,2)  =  0    (4-4) 

We  remark  that  no  new  identities  are  obtained  by  considering 
systems  consisting  of  more  than  three  particles.   Actually,  we 
get  no  nev;  information  even  when  we  consider  a  system  of  three 
particles  if  we  are  in  case  (a),  or  if  we  are  in  case  (b)  and 
assum.e  that  the  intermolecular  forces  are  central. 
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One  can  verify  the  following  formula  of  Irving  and  Kirk- 
wood  [15]  by  direct  differentiation: 

If  X,,  Xp  are  vectors  and  K  is  a  vector  function  defined 


by 


where  *  is  some  smooth  function  and  where  the  Xp-integration 
extends  over  the  whole  space,  then 

div  K  =  ^.  K  =  |y[<D(x^,X2)  -t.(ir2,ir^)]dX2  (4-6) 

Conversely,  the  right  hand  side  of  (^-6)  is  the  divergence  of 
K  (as  defined  by  (4-5))  plus  an  arbitrary  vector  function  whose 
divergence  vanishes. 

We  shall  make  repeated  use  of  this  formula  to  identity 
certain  terms  as  divergences  and  the  corresponding  flow 
vectors  with  vectors  of  the  form  (4-5) .   To  verify  that  the 
latter  identification  is  reasonable,  we  have  to  'calculate  the 
flow  of  e  across  the  surface  of  a  domain  D.   In  order  to  do 
this,  it  is  necessary  to  localize  the  flow  of  e.   For  the  part 
of  e  which  is  carried  by  the  molecules,  this  has  already  been 
done  by  the  localization  of  e.   For  the  part  due  to  inter- 
molecular  forces,  it  is  natural  in  case  (a)  to  localize  the 
flow  between  two  molecules  along  their  line  of  centers;  in 
case  (b),  it  is  not  so  natural  (since  the  force  of  molecule  1 
on  molecule  2  can  transfer  e  to  molecule  3),  but  justification 
can  be  given  (see  Ref.  [11]),  and  it  is  assumed.   Without 
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actually  carrying  out  the  computation  (see  [15]  for  an  example), 
we  state  the  result  that  the  flow  vector  differs  from  a  vector 
of  the  form  (^-5)  by  one  which  is  small  and  localized  and 
whose  divergence  vanishes. 
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5.    Conservation  Equations  of  Fluid  Dynamics 

If  we  multiply  (2-1)  by  Ne,  and  integrate  with  respect 

4.    dP 

to  -^^,    we  get 


de,    N    __         5e,      5e        Be        be 


-  (N-l){Fl2-^+°12-  35^)  =^  '  (5-1) 

where  we  have  assumed  that  f  vanishes  (or  is  periodic)  at  the 
limits  of  integration  of  all  the  variables. 

Similarly.,  if  we  multiply  (2-1)  by  Ne-,2  and  integrate  with 

dP 
respect  to  -r—  ,  we  get 


^e;L2^  a    ,f=      ,       r]        ^12^     ^^12^,    ^^12^     ^^12 
-?t-+^  •  (^1^12)  -  t^l*  -S3^+Pl'  -^^^2'   -^  +  P2'  -^ 


+  (^11-^V-  4|f -^('^11^12)-  4i7  +  (F22-*-W'  ^ 


Bsp  ^^10  ^ip       ^ip 


+  G23-^)=0  (5-2) 

Note  that  the  density  e  of  e  in  3-space  is  given  by 


e  = 


=  e^     in  case  (a) 


and 


e  =  e,  +-i — §—^£-,2     ^^^   <^^se  (b). 
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Also  let 


and  let  u,  be  defined  by 

dP. 


Pl^l  =  ^N  J   hh   d5^ 


Let 


^1  =  ^1  -  ^1 


u,  is  the  macroscopic  velocity  of  the  fluid;  c-,  is 
called  the  peculiar  velocity  of  particle  1. 

Next  v;e  show  that,  in  the  absence  of  external  forces, 
e  obeys  an  equation  of  the  form 

||+3l^-{eu,)+3|-.Q  =  0  ,  (5-3) 

where  Q  is  some  vector  (called  the  flow  vector). 

Case  (a) :  If  the  system  is  isolated  and  if  we  use  (4-1), 

equation  (5-1)  becomes 

If  +^  •  («"i)  +gl^  •  (^i)  -  ("-!) (^)  =  °  (5-4) 

But  (4-2)  gives 

(1,2)  --•|[(1,2)  -(2,1)]  . 
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Hence 


||+^-(eu^)+^-(^,)-^[(i:2)-(2,l)]  =0       (5-5) 


By  (-^-6),  this  is  in  the  form  (5-3),  with 


^a  =  cp^-N(N-l)K^  ,  (5-6) 


where  K  is  given  by  (^-5)  in  terms  of 

cl 


dXg 


Case  (b) :  Let 


*  -     (N-1) 
^12    ^1    2    ^12 


Then 


■X- 

e  =  e 


12 


Again,  we  assume  that  the  system  is  isolated.   Multiplying 
(5-2)  by  ^^^-"-^  and  adding  it  to  (5-^),  we  get 


M+^*(^^l)+3|--(^l^l2)  -(N-1)  (172)  -^^^-^y^-^U(l,"2;i)  +(27173) 


(N-1), ^^12  ,_      ,  ^^12  ^   ,  ^^12  „   ,ifl2^N   (N-l)r/r=K=^N  ,  /^^^^^ 
-    2^(-35^'^l  ■'-^'^2   -"-^  Pi  "'"■^•^2)  -^-^[(1.2,2)  +(2,l,r) 

=  0  .  (5-7) 
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Now,  eq.  {^-^) ,    together  with  the  symmetry  of  f,    gives 


(1,2,3)  +(2,1,3)  -  ^[(1,2,3)  +(2,1,3)  -  (2,3,1)  -  (3,2,1)] 


=  ^[(1,3,2)  +(3,1,2)  -  (2,3,1)  -  (3,2,1)] 
Use  of  this  and  eq.  (4-3)  reduces  (5-7)  to 

-  (^^-i)^^-^)[ (17372)  +(37i72)  -  (27371)  -  (37271)  ]  =  0    (5-8) 

Eqn.  (5-8)  is  in  the  form  (5-3),  with 


(t> 


Q^   =   c^.;^  -  N(H-1)K^  .mzimtSl  K^    ,  (5.9) 

where   K     is   the   same   as   before   and  K,    is   given  by   ('4--5),    with 

f^be   ^  ^e..^  ^^-iQ  ^^1-^  dP^dPpdP^ 

b  =  j'(-^-^12  +^-^12  +-3f^-''32  ■'^•°32'^123      d;^4  '^"'°' 

The  equations  of  conservation  of  mass,  momentum,  angular  momen- 
tiAm,  and  energy  are  special  cases  of  eqns.  (5-5)  and  (5-8). 
Since  they  will  be  needed  later,  we  shall  write  them  dovm 
explicitly. 

We  adopt  the  convention  that  if  a  quantity  has  subscript  i, 
it  is  evaluated  at  x.;  if  there  is  only  one  space  coordinate 
involved,  we  will  drop  the  subscripts. 
Mass:     e,  =  m 
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Eq.  (5-5)  becomes 

||.|,.(pu)=0  (5-11) 


Momentum:       £-.  =  m^^  . 


Eq.  (5-5)  becomes 


|.(p,i)  ,_ij(p.V)  .ig  =  0  ,  (5-ia) 


where  p'^^"  is  the  stress  tensor  and  is  given  by 


with 


P   ~  Pkin   Ppot 


i  J 
=  mc  C^ 


pij   =  -N(N-1)kJ^* 
^pot     ^    ''  a 


and 

^  «   i     dP^dPg 

*a(^l'^2)  "/''  ^12^12  dX3_dX2  • 

If  the  intermolecular  forces  are  central, 

^2  =  (x2-"4)^*^'''2'  ""l'^  ^ 
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(5-13) 


*  *  ,        ,      -I  •^^Tp 

where  F  is  some  function.   (Usually,  F  (  ix^-x.,  |)  =  —   ,   , 


where  r  =  |xp-x, | . ) 
Let 


y^  =  Xg  -  x^ 


and 

dP^dP, 


r      1 
n-L2^^1'^2^  "7^12  d3^ 


2 


dXp 


Then 


Angular  Momentum 

Since,  when  i,j,k  are  cyclic,  the  i-th  component  of  the 
angular  momentum  depends  on  the  j-th  and  k-th  components  of 
position  and  velocity,  we  will,  depending  on  convenience, 
denote  it  with  the  superscript  i  or  jk. 

eI^  =  m(xJeJ-xJe^)  +M.f  , 

1k 
where  ll^   is  the  angular  momentum  of  molecule  1  due  to  motion 

relative  to  its  center  of  mass. 

We  have 


^1^1    ^l^lkin  +^1^1   ^iPikin 
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and 


-^'^12  =  "^'^la  • 


Hence   (5-5)    takes   the   form 

-^+A(^iM)  -^A  e^sf -N(N-l)  -V  =  0  ,        (5-15) 


ax^    ^    ^     bx^    ^  ■"  bxi 


.ijk 


where  k^^   is  of  the  form  (4-5),  with 

^M- J^         dP^  dP. 


*'  =/tm(xJp^2  -x^pJg)  t-^-a^^^Jf^g  ^3^  .         (5-16) 


12    '  '  ' 

Energy        e-^  =  ^i-^  +  e^ 


e  =f  +  (^-1)  f   , 


t  t  t 

where  e,    is  the  kinetic  energy  of  molecule  1  due  to  motion 
relative  to  its  center  of  mass.   If  the  molecules  are  rigid, 
such  motion  reduces  to  rotation.   In  that  case,  if  (>:>     denotes 
the  angular  velocity  of  a  molecule,  its  internal  angular 
momentum  u  is  given  by 


/  =  I^j'coJ  , 


where  I  "^  is  the  inertia  tensor  of  the  molecule.   Its  internal 
kinetic  energy  is  -^1   '^cu  cd^  .   In  the  general  case 
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t  I  I 


1   2  1  ,  ^Ik   k 


and  eqn.  (5-8)  becomes 


ax^;        ^x^       ^      ax'3_ 


with 


q^  -  2  mc-^c^  +  ^ikin^l  +^1   ^1  ^   2    ^12^1  ' 


and 


Be,'"         dP^dPg 


*a=/("^VFl2+-5i-'^12)f 


12   dp^  "^12^^  12  dx^dx^ 


<J>,  Is  given  by  eqn.  (5-10)  . 
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6.   Conservation  Equations  for  Correlations 

\<le   first  derive  the  equation  for  the  joint  density  of  a 
pair  of  Integrals  e',  e".   Recall  that  there  are  three  cases: 

(1)  Both  integrals  are  of  type  (a). 

(2)  One  integral  is  of  type  (a),  and  the  ether,   of  type  (b), 

(3)  Both  integrals  are  of  type  (b) . 
Recall  also  that  if 

^12  "  ^1^2  ^^  ^^^^    ^^'^' 


^12  "  ^1^2  +^^2^^£{e23  ^^  ^^^^    ^^^' 


and 


^12  =  4^2  +-^-^(  4^23-^  4^{3>^  (N-2)^(N-3)  ^  .  ^^n^    -^  ^^^^  (3^^ 

then,  in  each  case,  Q, p  is  the  joint  density  (in  3-space)  of 
e'  and  e"  when  x^  and  Xp  are  far  apart. 

If  one  integrates  eqn.  (2-1)  — the  Liouvillfe  equation  — 
with  respect  to  all  but  r  particles,  one  gets  vjhat  is  called 
the  r-th  equation  of  the  BBGKY  hierarchy: 

Sf.o       r  ^fno         ^f  ^  r   3f,o 

^- •  ■  =" -a-^- •  • -•  t, -^^- •  • -•  V -.1=  (^- •  • -^^ -1 3 

1=1    1         -1        ijj=l    1      ■' 

±^3 


+  _^-..r.G..)+(N-r)^/(~^"-^    'F, 


12.  .  .r  r-fl. 
^i 


+  _^---'  ^-^.Q.  ,.)dP  ^,  =  0  (6-1) 

op.         ir+l-*   r+1  ^    ' 
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(It  is  assumed  that  the  system  is  Isolated.) 

Case  1:    Multiplying  the  second  BBGKT  equation  by  N(N-1)Q  „ 

dP,dP2  ^"^ 

and  integrating  with  respect  to  -^ — -^ — ,    we  get 

3Q-|  p    ^  -^    ^i     ;ii    

^t--^^*(^A2)  -^^*(^A2)  +^-(^2^12)  -^^-(^2^12) 


^z\  ae-       ^s^       T^  "67^  . 

-  -53?7-^1^2-^'Pl^2-^-^24-^'P24-  (^12  ^-^^12' 3^) 


2 


-  (^21^+^21'^)4-(^"2)[(F^3.^  +  G^3.^)s^ 


If  we  use  eqn.  (^-1)  and  assume  that  |xp-x-  |  >  a,  eqn.  (6-2) 
becomes 


-  (''-2)[(''i3-^+V357'^2-^(''23-^-^°23-bJ)4]  =  °  (6-3) 

The  set  over  which  the  assumption  that  |xp  -x.,  |  ■>  o   fails 
is  a  strip  of  width  2a  in  six-dimensional  space.   But  a  is  a 
small  parameter,  and  ultimately  v/e  will  take  a  limit  in  which 
0  -*  0.   It  can  be  shown  by  considering  the  initial  value  problem 
that  if  the  terms  dropped  are  small  initially,  they  will  remain 
small  over  a  finite,  macroscopic  amount  of  time.   Intuitively, 
this  means  that  if  tv/o  specific  particles  are  far  apart  initially, 
they  have  a  very  small  probability  of  coming  close  to  each 
other  within  such  an  amount  of  time. 
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Using  (^-2)  and  the  formula  of  Irving  and  Kirkwood,  xve 
can  write  (6-3)  as  a  conservation  equation: 

-^+^•(^1^12)  +33^'(^2"^12)  +B77-(^A2)  -^'(^2^12) 


-  N( 


(N-1)(N-2)(X.JC(1)+^.J^(2))  ^0 


where 


j/(l)  =/4'k(1)(x^,P2) 


dP. 
dZ 


with 


*(^^x^,X3)  =/(^. 


dP-i^dP^ 


13  '  6p7  ^13^^123  dx,dx. 


and 


,(2)  =  /;.k(2)(p  X  )  ^ 


with 


dPgdP^ 


(6-4) 


(6-5) 


2  "23^  123  dXgdx^ 


Case  2:    On  multiplying  the  third  BBGKY  equation  by  N(N-l)Q^p 

dP.dPgdP-. 
and  integrating  with  respect  to   ,  ^ — —  ,    we  get 
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^^•(-A2)-^•(^l2)-(^l^-Pl5J^t^^^^^23^ 


Ss' 


-  (N-3)(F24'^+^24'-3^)4 


-  ^ 2 -^(^24--^+^24--^+V-^+^3^  ^1  ~ 

We  observe  that  the  bar  denotes  integration  (with  respect  to 
everything  except  x  and  x^)    and  that,  under  the  integral  sign, 
the  diminy  indices  can  be  interchanged  at  will.   For  example. 


^24'  ^^1  -  ^23  "^^1 


If  this  property  and  the  identities  (4-1) -(4-3)  are  used  and 
if  |xp-x,  I  >  2a  then  eqn.  (6-6)  becomes 
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%+^-("l"l2)  +S;-("2«i2)  '■4l<'^12)   ^4^<-^12) 


dsi,  dsX 


5i^         5i;!         ^ 


(N-2)(N-3).^     1  4-n   Z-i  V  .^-r  .-J^W" 


-  ^^!=2KNz2i(,^,.!|a.a^,.!g2..3,.!|i.03,l||rei  =  0  (6-7)  ^ 

Use  of  {^-^)    -    (^-6)  enables  us  to  write  (6-7)  as  a  conservation 
equation: 


Sq 


12  ,  B  ,     -^ 


3t^+3f^-(Vl2)  +^'(^2^12)  +^'(^A2)  -^^-(^2^12) 
-  N(N-1)(N-2)(^.}C(1)  +^-3^(^b 


N(N-l)(N-2)  ^  .y(3)   N(N-lMN-2UN-3)  ^^  :y(6)  ^^   .g_Qs 
2      dx,  -^  3        <3xo  ^    "^ 


wherej(^^  ^  andJC,^  ^  are  the  same  as  in  (6-5)  and 


y3)  =y;^^K(3)(x^,P2,P3) 


dPgdP^ 
dXg 


with 


Se,'      Se,'  dP,  dP^^ 


'^^^  =A^*^14  -^■^•^14)^1234  d3^ 


dx^  ' 
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JC(^)-/eiK(6)(P^,.,)^ 


with 

Case  3:    On  multiplying  the  fourth  BBGKY  equation  by  K(N-1)Q^2 
and  integrating  with  respect  to  dx~dx ~  '  ^®  ^^* 

^+^-(-A2)  +^-(^12)  +^'(-2"^12)  ^^•(^M 


-(^l^-Pl5^)[^2-^¥^^23^ 

-  2-^(^1  -^^Pl^i^)'- ^2  +   2   ^24  J 
-(^2^-P2S|^t^{^^43^ 

-  -2^^^2'-^+P2'-^^  ^^1+2-^^13^ 

-  ^-2   ^^3'-^+P3'~^^^^2+   2-^^24^ 

-  -^-^(^4•-33^^-P4'^^)^^l+  2  ^13^ 

1=1   j=l    ^   ^x     ""    1  1         1 


[,.,^  ,iN^(e.e^3  .e^s'3)  .iii^^I^s-.s:',,  1 


4     ^13  24^ 
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-  td'si  ■'''32  +F3i,+(«-'*)f35)-5|^  +  («3i  +832  +«34  +(''-'*)°35'-^' 


-IC^l  +''42  +-"43  +(N-t)P45)'3%  +  ('54l  -^"ita  +°43  +(''-'*)°it5)'3%' 


[("-^V"-3)s'3B^,]  =0  (6-9) 

If,  as  before,  we  interchange  the  dummy  indices  under  the 
integral  sign;  use  identities  (^-1)  and  (^-3),  and  assume  that 
|xo-x,  I  >  30,    eqn.  (6-9)  reduces  to 


12  ,  ^   /   7T   \  ,   S   /„  -7T 


3T^-^Bf^-(-A2)  -^^%-(^2M  -^^'(^1^12)  -^^'(^2^12) 


-  2-^(^3  ^  +^13  ^  "  ^31  ^  ^31  ^>  ^2 

(N-2)  /p   ..   2^   „   .2   p   ._3_p   .  ,  i\pt 

2    ^^23  a^  ^^23  ^   ^32  ^    "32  ^^^1 

(N-2)(N-3)^p  .  ^23,p      23  .p   ._!23,p   ._!23n., 

-  2 ^(^24^T^+^24'-3iJ+^3^  -^+^-3^'^^^^l 


(N-2)(N-3).p   ^^13  ,r.    ^^13  .p  /^^Pip   /^13n.... 

-  ^    V    ^^  ^24-  ^  +  C^24'  3^  -  F^2'  ^  -  ^42-  ^)  43 

l^zmit3lf^       ^'14-r   ^'1  p   ^^3  r       ^'3x  .. 

-  ^ T-^(^i3  0^  -^^13  3^  ■  ^31*  ^  "  ''31  ^^  "24 


^e.U      ^e,'       Se'o      ^e-!. 


(N-2)(N-3)(N-4),„   ^"13.  p     13  ,p     13,  p     13Np.. 


-3^- 


=  0  .        '  (6-10) 

Use  of  Identity  (4-4)  and  of  the  Irving-Kirkwood  formula  enables 
us  to  write  (6-10)  as  a  conservation  eqiaation: 


#  -^•(^A2)  -^^•(^2Q"l2)  +^'("^12)  -^'(^12) 
-  N(N-1)(N-2)(^.}((1)+^,1C(2)) 

_  N(N-l)(N-2)(N-3).  5   ^^(3)  +^/^(4)) 
_  N(N-l)(N-2)(N-3)/  B  -u(5)  +J_v(6)^ 

3      ^^-K   ^^J^   ^ 
_  N(N-l)(N-2)jN:,3l(iMa(J-.3c(7)  -,^.X(8)^  ^  ,  ^      (g_,,) 

whereiC  ^"""^^Jn^^S-K^"^  ,  andj<^^  are  the  same  as  before,  and 

dP.  dP^ 


T^(^)  =  A'  K(^)rp    X    p  )  --1- 


with 


"^^  ^  =  j^^*^24  +^'^-24)^1234  dZ:d5^  ' 


3;  (5)  =Js^k(5)(x^,P2) 


dP^ 


dXg  ' 


^^;ith 
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dXg 


with 


,  .    >«ae'       ^e'       ae'       Be'  dP  dP  d] 


and 

dP,  dP, 


3C     J  13    i^i^^2'^3^   dx^  ' 


with 


*    =  j(-^'^25  -^^^'^25  ■'-^•^45  ■^-^•^45)^12345  dx^dx^  • 

Before  we  can  v/rite  dovm  the  conservation  equations  for 
correlations,  we  have  to  introduce  some  notation: 


^12  ~  ^12  "  ^1^2  ' 


^123  =  ^123-^2^3-^^23^^  ' 


and 


^1234  =  ^1234  -  hV3^i^  -  t^l^234^''  -  ^^1^2^34^^  -  «t^i2^34^^  ' 


where  a  is  some  number  and  where  the  square  brackets  indicate 
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a  symmetric  stun;  the  number  appearing  above  each  sum  indicates 
the  number  of  terms  in  it.   For  example. 


tV23^^  =  V23+Vl3"'^3^l2  ' 


The  above  formulas  serve  as  definitions  of  ■'/'-12' ^123'  ®"^<^*   ^'^ 
analogy  with  the  definition  of  the  r-particle  distribution 
function,  ^, ^     is  called  the  r-particle  correlation. 

The  notation  introduced  above  makes  it  possible  to  express 
the  r-particle  correlation  function  in  terms  of  the  r-particle 
distribution  function  and  lower  order  correlations  or,  equiva- 
lently,  lower  order  distribution  functions.   Thus  if  r-particle 
correlations  are  ignored,  the  approximate  or  truncated  r-par- 
ticle distribution  function  f, ^   „  can  be  expressed  solely  in 

X  ^  •  •  •  -L 

terms  of  lov;er  order  functions.   For  example, 

4\23  =  fif2^3^^V23^^  --^   -2f^f2f3+[f^f23]3. 
The  significance  of  a  is  that  specific  values  of  a  give  rise  to 
specific  ways  of  defining  correlations.   For  example,  a=l  gives 
the  well-known  Mayer  cluster  scheme.   For  more  details,  see  Sec.  7. 

If  e',  e"  are  any  two  integrals,  we  denote  the  correlation 
between  them  by  d^  ^  ^'^'^   its  convective  flow  by  a  ^  and  b   . 
Specifically,  in  Case  (1) 


dP,dP2 


12  =  ^A2  -  Vl   ^  =  N(N-l)yc^eje^V'i2  d^dx 


^2 
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b-^2        ^2^12 


~vr 


^r 


dP^dP^ 


e{   c^e^   =   K(N-l)Jc2BiB2-v!'^2   d^^dx^   ' 


in  Case    (2) 


,      _  TT      ^ "TT _  (N-2)~ rrr 

°12   ~   ^12"  ^1^2  2      ^1^23 


dP.dP 


=  N(N-l)/e's^^^2   d^^ 

N(N-l)(N-2)  r  ,    „    ,,,  .   ,^      ^^   ^      )    ""1""2" 

+  2  ^^1^23^^123      2^13      3   12^        dx-^dx 


dP^dPgdP^ 


dP.dP. 


■-12 


=   N(N-l)jc^£i£^V'-L2   d3^ 


dx. 


+ 


dP^dP^dP^ 


2  7/^1^1^23^^123^^2^13^^3^12^        dx^dx 


dP.dP. 


^2    =  N(N-l)yc2£-[s^^,2   d5^ 


dx. 


dP^dPodP 


.    N(N-l)(N-2)  r„    ^,     n    ,,  ,^   ,,,       ,^   ,,.      N         1      2 

+  -i y yc2e{e23^^l23^^2^13"^Vl2''        dx^dx 


3    . 


and   in  Case    (3) 


—  '~{     (N-2) 
-  pj  p^  — i — = — '- 


^12   ~   ^12  "  ^1^2  2 


■f  7^e~  +7^^'  )  -  ^N-2UN-3)—  ^ 
^^1^23^2^13^  ^+         ^^13^24 


dP.dP. 


=  NCN-DjTeje^'/^^g   d3^ 


+ 


dP-j^dPgdP^ 


N(N-l)(N-2)  r  .     „    ,  .  .„    ,       ,^   ^,^      N         1      2 

-^^ 2 :/^1^23^^123"^^2^13'^Vl2''         dx^dx 


dP^dP2dP^ 


+  2  J^2^13^^123      1^23      3^12''        dx^dx. 
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dP  dP2dP^dP;^ 
dx,dXp    ' 


vihe  re 


tV2'^34^^'  =  [fiV34l^- V3'^24-^2Vi3 


and 


[^12^34^^'  =  "f '^12^34]^  "^13^24  • 


a,  p  and  b  p  are  defined  in  the  obvious  v/ay. 

We  are  now  ready  to  write  down  the  equations  for  the 
correlations . 
Case  (1):   Using  eqns.  (5-5)  and  (6-4),  we  get 

^^^l^-tVia)  +J;-(Vi2)  +^1^-12  ^^-^2  +^iK<'*' 


where 


•j/(l)'=  -N(N-l)(N-2)j^^^^^  +i2N(N-l)K^^  , 


and 


j/(2)'  =  -N(N-l)(N-2)j/(2^  +epi(N-l)K2^ 
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(The  K's  occuring  in  eqns.  (6-12)-  (6-l4)  are  those  occuring 

in  eqns.  (5-6)  and  (5-9)-) 

Case  (2):   Using  eqns.  (5-5),  (5-8)  and  (6-8),  we  get 

-  4-*''' '-#^   =  0  (6-13) 

(2)  ' 
,  .1/(1)'    and  K'^    '      are   the   same   as   before,    and 

wherej^^    '  J\ 

.A3)'    _        N(N-l)fN-2)(N-3)-J3)    ,N(N-l)(N-2)         =7? 
Jn  ~  ~  2  J^        ^  2  la^23   ' 


and 


V(6)'    _        N(N-l)(H-2)(N-3)..(6)       N(N-l)(N-2)         rr" 
J\  ~  3  -^3  2b  1 

Case  (3):   Using  eqns.  (5-8)  and  (6-11),  v;e  get 


^+^-(Vl2)  +3lj'(Vl2)  -*-^'^l2+3l^*^2 

(6-14) 

,3l_.(3,(l)^3,(3)^J,(5)^.^<7)•^,^.(^(2).,^^(4).^^J6).,^^^^ 

flV        f2V        f3V  f6^' 

where  J(^    ^    ,X  >  J<.  y    andjc^^    >      are   the   same   as   before,    and 

^(4)'    _        N(N-l)(N-2)(N-3)..(4)      N(N-l)(N-2)         =t 
Jv  -  2  ,A  2  ^2a^l3   ' 

(5)'    _        N(N-l)(N-2)(N-3)^(5)    ,  rn   N(N-l)(N-2) 
X  ~  ~  3  ^^         '  ^2  3  ^Ib   ' 
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t/(7)  •  _  _  N(N-l)(N-2)(N-3)(N-4).(7)  +Ilh3lr^  H(N-l)(rJ-2) 


and 


(8)'  _    N(N-l)(N-2)(N-3)(N-4)  (8)  ,  (N-3)r^  N(N-l)  (N-2) 
_K     6        ^K   +  2  ^13      3     2b 

Finally,  we  shall  write  down  the  special  cases  of  eqns. 
(6-12) -(6-l4),  when  the  integrals  are  mass,  momentum,  angular 
momentiim,  or  energy. 
(1)  Mass  -  Mass 


e'  =  m,    Ep  =  m  . 

2  r         dP  dP 
%2  =  N(^-1)^>12  dF^-dif 

2  r  i      cip  dPp 

^12  -   N(N-l)mJc^V'-L2  d5^i^ 

.  p  ^  .     dP-,dPp 

^12  =N(N-l)myc^^^2^^  . 

-1^  "^  and -1/^ ^  """  vanish.   (6-12)  becomes 


Sx^-  -^  "■"       dx^     '  -''^    3x^   Sx^ 


(2)   Mass  -  Momentum 


E^   =   m,    £2*^  =  m?,^ 
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P  >»      dP  dP„ 
dP.dP. 


^2  =  ^^i^-i>^foi4h2  d5^ 


dXg 


b^^  is  defined  similarly. 


J/(1)'1J  =0  , 


and 


J;(2)'iO  =  -N(N-l)(N-2)j^(^)=^j'  +p^N(N-l)K^^  . 

with 

/o\  •     r  -;      dP^dP^ 
4.     -  mjF23f-L23  dxgdx^  ' 

Let  y'^  =  x^-x^.   Then,  if  the  intermolecular  forces  are 
central,  F^^  depends  only  on  r  =  |y|: 


F^3  =  y^'F*(r) 


Let 


dPgdP^    , 


^123  dXgdx^ 


Then 
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'123  ' 


and 


where  we  have  suppressed  the  dependence  of  Z-^2'\   °^  "^i' 
Recalling  that 

4i  =  f /yVF''(r)j^  n23(x2-Ay,X2+(l-A)y)dMy  , 
we  obtain 


dP, 
» 


p  "I  d  P 

^(2)  .ij  ^  .  ""("-lH^-g)/yVF*(r)/^(.A^g3+f,»,3tf3»,g)dXciy  ^ 

where  the  integrand  of  the  inner  integral  has  already  been 

dPgdPo 
integrated  with  respect  to  ^^  ^^r^   and  is  evaluated  at 

(x2-Ay,X2+(l-A)y)  , 

Eq.  (6-12)  becomes 

(3)   Mass  -  Angular  Momentum 


e'    =  m,    e"J^  =  m(xgep  -  x^d)  +u-^^ 


p  —      ilH^  ApC,p  Apl-ipy       1    M.p 

dP.  dP^  /»  .,  dP^dP, 


2 


dXg 
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T  "1  Ic         *i  "1  Ic 

a,x  and  b^M  ^^s  defined  in  the  usual  way. 


j/(l)'ijk  ^  Q 


and 


j^(2)  'ijk  ^  _N(N-l)(N-2)j/(^)^'^"^ +N(N-1)P^K^J 


a 


Eq.  (6-12)  becomes 

(4)   Momentum  -  Momentum 


ej^^  =  m^l  ,    £2  =  m^2 


dP.  dP^ 


dg  =  N(N-l)m2y?^?^f^2  ^ 


dXg 


If  the  intermolecular  forces  are  central,  v;e  proceed  in 
the  same  way  as  in  the  mass  -  momentum  case  and  obtain 

P  1  dP 

^(1)  'ijk  ^  _  m  N(N-l)  (N-2)|yV^k^.-(,)|J,^^^,f^^^^,f^,^^),,dy  ^  , 

(y^  =  x^-xj,    r  =  |y|) 


and 
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;  /ill  I     I X 

(y    =  ^3-  ^2'       ^  "  'yl^ 

Eq.    (6-12)    takes   the   form 

"^^12    ,     5    .    i,jkx   _^    a    ,    i,jk.    ^"^^12      ,'^'^12      ,^K ^^X ,. 

~st~  — i^^i  12^    — j-(^2  12^    — r~  — T"^^ — I —  1 —  ^ 

Sx-.  dxo  by--.  bxp  ^x,  Sxp 


(6-18) 


(5)   Momentum  -  Angular  Momentum 


i      i    "  '  1 


"  '  1 
where  e-,  "^  is  the  momentum  of  molecule  1  due  to  motion 

II  I-,- 
relative  to  its  center  of  mass.   Clearly,  e-,  "^  =  0  if  the 


molecule  is  rigid, 


sf^=m(x^^2^-x|4^)+U^^ 


Defining  d^p  -.  a^^  ,  and  b^^   in  the  usual  way,  we  see  that 
-^  ax^^'l'l^  ^   ^^"2^12  )  %^i   -^   ^  ^    ^ 

+  ^=^1 ^ =  0  (6-19) 

Sxp 
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(6)   Angular  Momenti.xm  -  Angular  Momentum 


e£J^  =  m(xJ^^-x^ej)  +Lif 


e^        =   m(x2^2  "^2^2^   ^2  ' 
Defining  dg^"^,  a];^^"^"^,  and  bj^^-^'"^  as  usual,  we  get 

.  (2)'ijkM 
+  -2^ J =:  0  (6-20) 

(7)   Mass  -  Energy 

e-I  =  m,    e'^  =  m  -^  +  e^^  ^ 


dP,dP„         A       dP.dP, 


d,2  =  N(N-l)|-/e2\2  d3^+N(N-l)"y^2''^^2  d5^ 


2 


dX2 


,  N(N-l)(N-2)mr„  ,,  ,  „  ,   ,  .  .   s  dP^dPgdP^ 

+  -^ ^i j'^23^^123+Vl3'^Vl2^    dx-^dXg 

a-|p  and  b,p  are  defined  as  usual,  and 


yi)'i  =  o=j^(3)'i 


So,  eqn.  (6-13)  becomes 


(6-21) 


^t--^(Ul^2)-^(Vl2)-^  -^'^^    -^^    ' 


(8)  Momentum  -  Energy 


We  define  d^p^  ^tp'  ^^'^   ^t"?  according  to  the  formula  for  case 
(2).   Eqn.  (6-13)  becomes 


+  _l-.(^(2)  'ij  +^(6)  'ij>)  ^  0  (6_22) 

OXq 


(9)  Angular  Momentum  -  Energy 

e 


e 


2 


e^  =  m-^  +  e^" 


The  eqn.  for  the  correlation  is 


+  _^(  (2)'i0k^_^(6)'ijk^  ^Q  (g_23) 
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(10)  Energy  -  Energy 

.2 


1    2   '  ^1 


.2 

^2     2    ^2 


d, p,  a^pvand  b'o  are  defined  according  to  the  formulas  for 
case  (3).   Eqn.  (6-l4)  looks  like 

^"^12,  a,  1^   W-A-fu^d  )    ,^"•12^^12^  ^(..(l)'i+..(3)'i^p{5)'i 
+^(7)'i)^_A^(^(2)'i^_^(4).i^(6)'i^_^8)'i^  ^0     (g_24) 


.48- 


7.    Generalization  of  the  Boltzmann  Equation^- 

The  Boltzmann  equation,  which  describes  the  evolution  of 
f,  for  a  rarefied  gas,  is 

-^+^1'^  =  (N-1)  /(fi.fg.  -  t^r^)\^dc.^^^i^  (7-1) 

The  integral  on  the  right  hand  side  is  five- fold  and  is  called 
the  collision  integral.   Primes  denote  quantities  after  collision; 
the  colliding  molecules,  namely,  molecules  1  and  2,  are  assumed 
to  be  at  the  same  point  physically.  V-.^   is  the  relative  speed 
(before  and  after  collision)  and  dcu  g  is  an  element  of  area 
in  a  plane  perpendicular  to  the  relative  velocity  ^2  ~  ^1' 

There  exist  several  derivations  of  the  Boltzmann  equation. 
For  quick  reference,  we  shall  give  below  a  rough  sketch  of  one 
of  them,  which  is  due  to  Grad.   (For  a  complete  account  and  a 
critique,  see  ref.  [1].)   All  the  derivations  make  two  basic 
assumptions: 

(1)  Binary  Collisions:   The  only  significant  collisions 
are  binary  collisions,  contributions  due  to  higher  order  colli- 
sions being  negligible. 

(2)  Molecular  Chaos:   Two  molecules  which  are  approaching 
are  uncorrelated. 

In  addition,  a  limit  process  is  Involved  and  it  is  assumed 
that  f-,  doesn't  vary  much  over  the  dimensions  of  a  collision. 
Following  Grad,  we  start  v.'lth  the  Llouvllle  equation  for 


The  discussion  in  this  section  pertains  only  to  rarefied  gases. 
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an  isolated  system  with  an  intermolecular  potential  which  is 
central: 

M+l^B^-eiti!:  #-^ij  =  °  (7-2) 

1=1   1    1,  j=i  ^1  "^ 

Let  p  >  0  be  some  niimber,  and  let  ro   denote  complements. 
For  i  7^  r,  let 


and 


D^^^  =  [P  :  Ix  -X.  I  >  o,    all  i    ] 
r     '-  r   '  r   1 '    '  ^r' 


S^^)  =  [x  :  |x  -X.  I  =  a] 
r     "^  r   '  r   i '     ■' 


Let 


D  =  D^-'-^x  Di^-'-^x.  ..xD^-^'  , 


^l^^l)  =yfdP2dP3...dPj^  , 


and 


^12(^1' ^2^  ^/^^^S  dP^.-.dPj^ 


d(1)xD^-'^)x..xd(^1) 


f^  gives  the  probability  of  finding  molecule  1  without  any 
neighbors  vrithin  a  radius  of  a.      f^^,  on  the  other  hand, 
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gives  the  probability  of  finding  molecule  1  with  at  most  one 
neighbor  and  is  not  symmetric  in  molecules  1  and  2,  since  the 
domain  of  integration  depends  on  1  but  not  on  2.   Also 

_  /»         N-1   /« 

fl(Pl)  =  fi(Pi)  -  (N-l)j  f-L2dP2  +(    2    U      ^ISS'^VPa  -•••  .. 

+  (-1)^"-^  I    fdrgdr^-.-dFj^    ;  -  - 

%^)    xd(1)  X...X4I) 


and 


fl2(Pl'P2)  =  ^12(^1' ^2^  -(N-2)  /  f^^^dP^+... 

+{-lf-^f    fdP3dP^...dP^ 

^(1)  x^^l)  x...xS(l) 

If  Na--^  is  small,  f^  and  f^o  approximate  f  and  f-,p  respectively. 
In  order  to  derive  the  equation  obeyed  by  f , ,  we  will  need 
the  following  formulas: 

/  -T — •A(x,,Xp)dx2  =  -  o  A(x,,Xp),d2   (Gauss's  theorem)    (7-3) 
Ixg-x^l  >  a  Ixp-x^l  =  a 


and 
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J   ^.A(x^,X2)dX2  =  ^'J  k{x^,x^)dx^  -ylAix.^,x^)'dZ  (7-4) 

Ixp-x-,  I  >  a  Ixg-x  I  >  a     |xp-x-,  |  =  a 


Now,  integration  of  the  Liouville  equation  over  the  domain  D 
gives 


-^ +y  I     -^'i-    dPo...dP„+^      /  -J^.F.  . 


dP^.  .  .dPj^  =  0 


Repeated  use   of    (7-'^)    gives 

f^'h   dP2"-dPN   =^-(Vl)    +(^-1)  //     Vl2(PrP2)-df  dt 

^2 

Similarly,  for  i  >_  2,  use  of  (7-3)  gives 


I^-k  ^P2---dPN  =  -//  Vi2(PrPi)-d^i  d^i 


"  ■//  ^12(^1' ^2)  ^2 '^^2  ^^2 


s(i) 

If,  as  usual,  we  assume  that  the  distribution  functions  vanish 
at  the  limits  of  integration. 


=  (N-1)  /--^'F^2  dPl^P2  • 


D 
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So  we  finally  get 


+  (N-1)  r   -^.F^2^P2  =  0  .  (7-5) 

^2 

If  we  now  choose  a  to  be  the  range  of  the  Intermolecular  force, 
the  second  integral  in  (7-5)  vanishes,  since  I^-io"'-'  when 
|xp-x-,  I  >  a.   So 

-^  + ^135^7  + (N-1)//   fi2(PrP2)(V^2)'^^2  ^^2  =  0       (7-6) 

s(i) 

^2 

This  equation  already  looks  a  bit  like  the  Boltzmann  equation. 
The  resemblance  can  be  made  closer  by  means  of  a  coordinate 
transformation.   Consider  the  diametral  plane  of  SA  '    v;hich 
is  perpendicular  to  the  relative  velocity  in'^i  '      '^^^   points 
in  this  plane  be  represented  by  polar  coordinates  (r, e)  where 
0  <  r  <  a  and  0  <  e  <  27r.   The  area  element  dco   in  this  plane 
is  given  by 


dcD^P  ^  rdrde 


Consider  the  projection  of  points  on  Sa  ^    onto  this 
plane.   The  plane  is  covered  twice,  but  the  mapping  can  be 
made  one-to-one  by  distinguishing  the  tv;o  hemispheres  formed 
by  the  plane: 
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Spi   (^  -^  ).d2  >  0   (receding  particles) 


Sp:   (^2~^i)*'^^  "*  "^    (approaching  particles) 


If  (r, e)  is  a  point  in  the  plane,  let  V^   and  PZ  respectively 
denote  the  points  on  Sp  and  Sl  which  map  into  it. 
Then  equation  (7-6)  becomes 

-^^^1^=    (N-l)/[fj2(PrP2)  -f?2(Pl'P2)^^12<^^12^^2       C^'^) 

Equation  (7-7)  can  be  transformed  into  (7-1)  by  performing 
four  steps: 

(1)  Replace  the  arguments  of  f^^^i^i' ^p' ^"^    ^^   ^^®  values  they 
must  have,  say,  P',P',t'  at  the  start  of  the  binary  collision 
in  order  that  they  may  have  the  values  P-i »  Po  ^-'^  time  t.   (The 
start  of  a  collision  is  defined  by  [x-j-Xol  =  o.) 

(2)  Put 


fi2(PrP;)  =  fi(Pi)f2(Pi) 


and 


^12(^1' ^2'''^'^  =  fJ(Pj^,  t')f2(P^,  t')   (chaos  assumption) 


(3)   Replace  xjjX'jxZ  by  x,  and  ti  by  t    (slowly  varying  f) 
(^)   Replace  f^"  by  f^  and  f^   by  f^    (Na^  small). 

It  follows  that,  in  oi-der  to  show  that  the  Boltzmann 
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equation  describes  a  rarefied  gas,  one  has  to  verify  the  two 
basic  assumptions  made  In  its  derivation.   It  is  intuitively 
clear,  though  it  hasn't  been  rigorously  proved,  that  the 
assumption  of  binary  collisions  becomes  asymptotically  valid 
as  the  gas  gets  more  and  more  rarefied.   Concerning  the 
assumption  of  molecular  chaos,  however,  there  seems  to  be 
some  confusion.   As  pointed  out  earlier,  it  is  often  stated 
in  the  literature  that  the  chaos  assumption  is  also  related 
to  the  density  of  the  gas;  indeed,  that  it  is  violated  only 
when  the  density  increases.   But  we  know  from  the  results  of 
the  last  section  that,  even  for  a  rarefied  gas,  there  exist 
situations  —  such  as  turbulence  —  in  which  correlations 
are  large  and  persist  for  microscopically  long  times.   In 
order  to  t3ke  such  situations  into  account,  we  propose  the 
following  simple  and  natural  generalization  of  the  Boltzmann 
equation  (actually,  of  equation  (7-7),  before  making  the 
chaos  assumption) : 

4t  +  ^1  ^  =  (N-l)/(f^,2'  -^12)\2n2d^2  ' 

^f-io       ^flO       ^flQ        ^^19        ^^10 

(N-2)j'(f^g,2,  -  fi23)V23*J^23de3  +(N-2)y(f^,23,  -  f^23  '^13*'°13*^3  ' 

and  so  on.   Here,  as  in  the  Boltzmann  equation,  we  have  assumed 
that  binary  collisions  are  the  dominant  ones  and  that  the 
colliding  molecules  are  at  the  same  point  physically. 
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The  first  two  equations  of  the  sequence  (7-8)  have  been 
postulated  and  used,  together  with  truncation  of  f-ioo^  fo^ 
the  calculation  of  transport  coefficients  by  Devanathan  and 
Bhatnagar  [lb  ]  . 

The  sequence  (7-8)  can  be  made  plausible  in  a  manner 
similar  to  Grad's  derivation  of  the  Boltzmann  equation.   Note 
that  the  first  equation  is  what  the  Boltzmann  equation  would 
be  if  the  chaos  assumption  v;asn't  made.   To  obtain  the  second 
equation,  for  example,  let 

E^  =  D(.l)n  d(.2)    (r=3,4,..,N)  . 


Let 


fl2(PrP2)  =f^   dP3dP^...dP^  , 


E^  X  E^i  X  .  .xE^ 


and 


fl23(PrP2'P3)  =f^   dP^dP^---dPj^  • 


Ei|  X  Ef-  X  •  •  •  xE«T 


..oo 


f^P  is  symmetric  in  particles  1  and  2;  f-,po  is  invariant  under 
the  interchange  of  1  and  2  but  not  of  1  and  3  or  of  2  and  3- 
In  other  v/ords,  f^po  gives  the  probability  that  the  cluster 
consisting  of  molecules  1  and  2  has  at  most  one  neighbor.   Also 


If 

They  claim  that  the  equations  describe  a  moderately  dense  gas. 
¥e  believe,  however,  that  the  dense  gas  correction  is  incomplete. 
The  clearest  indication  of  this  is  the  fact  that  the  conservation 
equations  as  well  as  the  equilibrium  f-^p  they  (the  sequence  (7-8)) 
lead  to  are  those  of  a  dilute  gas. 
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N-2 


flg^Pl^Pg)  =  fi2(?r  ^2)  -  (N-2)J^f^23^P3  +(    )/ f  i234^P3dP4 


E-3  E^xEn 


Eo  X  El.  X  .  . xE,^ 


and 


f\^^{v^,v^,v^   =  f^23(PrP2'^3)  -(N-3)/fi234'^P4  +  --- 


•V    % 


+  {-l)^~^r  f   dP,,dP^.  .  .dP. 


Ei|  X  E(-  X  .  .  .  xE,y 

It  follov7s  that  if  N  —  0  and  a  —  0  such  that  Na^  -»  0  then  f^g 
and  f-,po  approach  the  full  distribution  functions  f,  ^  and  f-joQ 
respectively. 

Integration  of  the  Liouville  equation  over  E^  x  E^^  x  .  .  .xEj^ 
and  use  of  the  same  type  of  reasoning  as  in  Grad's  derivation 
(except  for  the  chaos  assumption)  gives  the  second  equation  of 
the  sequence  (7-8). 

A  function  (of  velocity)  is  called  a  collisional  invariant 
of  a  kinetic  equation  if,  on  multiplying  the  collision  term  of 
the  equation  by  the  function  and  integrating  with  respect  to 
the  velocity,  one  gets  zero.   The  word  "invariant"  is  proper, 
since  the  operation  referred  to  gives  the  average  change  (per 
unit  time)  of  the  function  due  to  collisions.   The  collisional 
invariants  of  the  first  generalized  Boltzmann  equation  (G.B.E.), 
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i.e.  of  the  first  equation  of  the  sequence  (7-8),  are  the  same 

i  2 

as  those  of  the  Boltzmann  equation  (B.E.),  viz.,  1,  ^,  (i=l,  2,  3) ,  ^-|  . 

The  proof  is  the  same,  too: 

Let 

I  =/(^1.2'"^12)^12'^^12'^^2 

=  /t^l2(^i'^2)  -fi2(^i'^2)^^12^''l2^^2  ' 


and  let  *  be  any  function  of  velocity. 
Let 

I^)  =/^(^i)(fi.2.  -fl2)^12^^12^^2^^1  • 

Interchanging  C-i  and  ^p,  and  therefore  ^'  and  il,    and  observing 
that  the  Jacobian  is  1,  doD^2=dcUp- ,  and  V,  p=Vp-|,  we  get 

\   =  \  =/*(^2)(fl.2.  -  fi2)^12^-12d^2<^^l 

Here,  we  have  used  the  symmetry  of  f^p  in  the  two  particles 
as  well  as  the  assumption  that  the  two  particles  are  at  the 
same  point.   Similarly,  interchanging  ^^  and  ^',  and  Co  and  C', 
we  get 

*1     *1     *2 
where 
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i<t,i  =/*(^i)(fii2' -^12)^2^^12^^^!  ' 


and  similarly  for  I.,  . 
So 


A-i 


\      =i/[*(ei)  +*(?2)  -'^(4)  -*(^2)^t^l-2'  -^12^\2^^12^M^1  • 

(7-9) 


Hence 


*1 


if  *  is  a  summational  invariant,  i.e.,  if 


M^-l)  +^(^2)  =  ^{i{)   +*(^^)  (7-10) 


By  virtue  of  the  laws  of  conservation  of  mass,  momentum, 

1   2 
and  energy  this  is  true  of  1, i^,  i^    .      By  a  theorem  of  Grad, 

any  continuous  function  *  satisfying  (7-10)  is  a  linear  combi- 

i   2 
nation  of  1,  ^. ,  C-,  • 

We  remark  that  the  reason  for  the  similarity  between  the 

above  proof  and  that  for  the  Boltzmann  collision  term  is  that 

f^2  has  the  same  symmetries  as  f-ifo*   Also,  the  same  type  of 

reasoning  as  above  shows  that  the  collisional  invariants  of 

the  second  G.B.E.  are  products  of  those  of  the  first: 
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wnere 


a,f3  =  0,1,2,3,4  , 


and 


¥e  have  not  been  able,  however,  to  show  that  all  the  collisional 
invariants  are  of  this  form. 

The  sequence  (7-8),  like  the  BBGKY  hierarchy,  has  the 

*  ( r)  fr+ll 

property  that  the  equation  for  f^  '    contains  f^    '  in  the 

collision  integrals.   In  order  to  solve  the  sequence,  therefore, 

one  should  terminate  it  at  some  point.   If,  for  example,  it  is 

desired  to  terminate  it  at  the  r-th  stage,  then  in  the  equation 

(t)         fr+ll 
for  f^  '',  f^    '    should  be  replaced  by  its  truncated  form 

/Vr+ll 

f^    ',    which  depends  only  on  the  lower  order  distribution 

functions  f^  ^,f^  \  . . . , f^^' .      Clearly,  the  truncation  cannot 

be  done  arbitrarily;  certain  compatibility  criteria  should  be 

met.   One  obvious  compatibility  condition  is 


/^12...r   dPi  =  f^"'=A2...r   ^^i    (i=l.  2,  .  .  .,r) 


CI 


A  second  condition  is  that  if  the  tmncated  r-th  G.B.E. 
is  called  the  r-th  B.E.,  one  should  have: 

Integration  of  the  r-th  B.E.  with  respect  to  any  particle 
yields  the  (r-l)-th  G.B.E.  02 


•X- 

When  the  subscripts  are  unimportaJit,  we  shall  denote  the  r- 
particle  distribution  function  by  '>(r). 
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The  first  condition  means  simply  that  the  truncated  function 

satisfies  the  same  normalization  condition  as  the  untruncated 

(r)      Afr") 
one;  in  other  words,  the  correlation  f  ^  -"  -  f  "^  ^    vanishes  when 

integrated  ^vith  respect  to  any  of  the  particles.   The  second 

condition  guarantees  that  trimcation  at  the  r-th  stage  doesn't 

imply  any  loss  of  information  at  the  lower  levels.   It  may  be 

worth  noting  that  several  recent  generalizations  of  the  Boltzmann 

equation  [5,17,18,19]  satisfy  not  C2y  but  the  more  restrictive 

condition  that  the  r-th  B.E.  gives  the  (r-l)-th  B.E.  (not  G.B.E.) 

on  integration. 

There  are  several  ways  of  truncating  the  distribution 
function  f^^^  subject  to  conditions  CI  and  C2,  but  we  shall 
discuss  only  three  of  them  in  detail.   If  we  impose  a  third 
condition,  the  three  truncations  turn  out  to  be  equivalent  in 
a  certain  sense. 

Consider  the  moments  about  the  mean  of  a  random  variable 

X. 

Let 


x=/. 


X  dp  , 


where  dp  is  some  probability  measure, 
For  any  positive  integer  n 


/(x-x)^dP  =  x^-nx^-1  x  +  (g)x^-2x2-..  +  (.i)^-l(n-l)~^ 


Making  the  association 
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f(l) , 


j{x-~)^   dp-*  i'^^'    (the  n-particle  correlation  function)  , 


and 


x^  -  f  (^)  , 


we  obtain  a  scheme  suggested  by  Grad: 


f  12  =  V2  +  ^2 


^123  =  V2^3''fV23^^+^123 


^1234  =  V2Vi^  +[fi^234^^  +^^1^2^34^^  +^1234 


and  so  on.   Here,  as  before,  the  square  brackets  indicate 
Slims  symmetric  in  the  subscripts,  and  the  number  appearing 
above  each  sum  indicates  the  niomber  of  terms  in  it. 

Similarly,  by  considering  the  cumulants  of  a  random 
variable,  one  can  obtain  Mayer's  scheme: 


^12  ~  ^1^2  "^^12 


^123  "  ^i"2^3  "^  f^  ^1^23  ^  ■*"^123 


^1234  =  ^i^2^3^4  +^^1^^234^^  +  [f^f^[i^^  f  +  [u-j^nU-^i^]^  +^^1234 
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and  so  on.   In  general 

^123... r  =  V2---^r  +  f  ^1^234... r^+tfi^^23^45.,.r]+[flM.234^56...r^ 
+  ...  +  [flf2^^34...r]+[V2^34^56...r^+••• 
+  [V2---V2^^r-l,r^  +  ^^12^3^.  .  .r  ^ 
+  [p.  „„[Xk(-    ]  +  all  distinct  products  of  the 

^t^+^^234...r  • 
Finally,  one  can  also  write 


^12  =  ^1^2  "^^2 


^123  =  ^l^2^3  ■*"  [^1^23^  "^^123 


^1234  =  fiV3^4  -"f  ^1^234  ^"^^f  ^12 V^^  ''^234 


and  so  on.   In  general 

^123. ..r  =  V2---^r+tV234...r^+ffl2V...r^+t^l23S...r^ 

+  ...  +  [fi234...r-2^r-l,r^  +\2...r  * 

The  correlation  functions  in  the  three  schemes  are  related 
to  one  another  as  follows: 
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^12  =  t^l2  =  \2 


^123  ""  ^123  ~  ^123 


L    (7-11) 


3  _ 


^1234  "  ^1234  +t(^i2^'-34l  =  ^234  +^t^i2^34^ 


10 


,10 


^12345  ~  ^12345  "^^^123^45^   -  ^12345  "^^^^123^5^ 


and  so  on. 

The  truncated  r-partlcle  distribution  function  is 
obtained  by  simply  setting  the  r-particle  correlation  equal 
to  zero.   For  example,  the  truncated  four-particle  distribution 
function  is  given  by 


^1234  =  ^i^2^3^4  +^^1^234]^  +[f^f^i^-^i^]^ 


(Grad's  scheme) 


^1234  "  fif2^3^4  "^f ^1^234^  +  [i'lf2'^32.]^  +[m--l2M-34]^  (Mayer's  scheme) 


A      ^  4 

■^1234  ~  "^l'^2'^3'^4  ■*"t^i'^234^  "*" '■^12^34  ■' 


(the  third  scheme 


(The  truncated  one-,  two-,  three-particle  distribution  fujictions 
are  the  same  in  the  three  schemes,  since  the  correlation  functions 
are  the  same.) 

The  three  truncations  can  be  combined  in  one  formula: 

^^1234  =  ^i^2^3^4  +[^234^^ +[^1^2^34  ^^'*"°^[  ^12^34^^  ' 
where  a=0  corresponds  to  Grad's  scheme,  a=l  to  Mayer's  scheme, 
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and  a=2  to  the  third  scheme. 

Similar  formulas  for  higher  order  truncations  can  be 
given,  but  they  are  more  complicated.   Suffice  it  to  say  that 
the  truncated  r-particle  distribution  function  is  given  by 

^12. ..r  =^l^   ^i+tV23...r]  ^^^  ^23^45.  . .  r  ^  +  ^^  fi^23^45%.  .  .  ^] 
+  ...  +7[fi^234^6...^]+...+  [fif2%5...r^ 
+  5[f^f2^34^5g^^^^]  +  ...  +  [f^f2...fr_2Vl,r^ 
+^[^n2^-^4   r-'  "^  ^^^   distinct  products  of  the  f   . 

Proper  choice  of  the  coefficients  a, P...  gives  Grad's 
scheme,  Mayer's  scheme,  or  the  third  scheme.   But  if  we  let 
them  be  arbitrary  constants  subject  to  CI,  we  get  what  may  be 
called  the  polynomial- type  truncation.* 

All  the  truncations  mentioned  above  satisfy  CI  and  C2. 
We  shall  verify  this  for  Grad's  scheme,  but  the  proofs  are 
similar  in  the  other  cases. 

Note  that 


/f^2  dPg  =  /f^fgdPg  =  f^ 


The  Kirkwood  superposition  hypothesis,  which  states  that 

^  f  -p 
^        -   "Ig  23  13 

123  -   fif2^^3"  ' 

and  its  generalization  [20]  constitute  a  different  type  of 
truncation.   In  the  limit  of  small  two-particlle  correlations, 
the  Kirkwood  hypothesis  coincides  with  the  three  types  discussed 
here.   No  such  relation  seems  to  exist  between  its  generalissition 
and  the  three  schemes  in  the  higher  orders. 
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I.e. , 


fh2^^2   =--  0 


By  symmetry. 


A2'^Pl  =  0 


Using  this,  one  can  show  that  f-,no   gives  zero  when  integrated 
with  respect  to  any  of  the  particles.   For  the  general  case, 
let  Tp^'^',    Tp^^',    ...  ,-^^^~'^'    have  this  property.   Then  if 


^12.. .r=,^,  fl+[V23...rl'  +  [V2*34...rl'2'  + 


+  tfif2-.-fr-2*r-l,r]*2^ 


is  integrated  with  respect  to  any  particle,  say  the  r-th  one, 
the  only  terms  in  the  symmetric  sums  which  give  non-zero 
contributions  are  those  in  which  the  subscript  r  is  not 
associated  with  the  correlation  functions.   One  gets 


2J  fl+*12...r-l+[V23...r-l''"'+tV2*3t...r-ll*''2  ^ 
+  [fif2---fr-3*r-2,r-ll''"2  '  , 


v/hich  equals   f-,^        ^-1 
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We  observe  that  the  r-th  G.B.E.,  '.vhen  it  is  integrated  with 
respect  to  a  particle,  gives  the  (r-l)-th  G.B.E.   To  shoxv  this, 
one  has  to  use  the  same  argument  as  is  employed  to  find  the 
collisional  invariants  and,  in  addition,  assume  either  that  the 
force  terms  on  the  left-hand  side  of  the  r-th  G.B.E.,  which 
contribute  on  integration  to  the  collision  terms  of  the  (r-l)-th 
G.B.E.,  can  be  put  (after  integration)  in  the  form  of  the 
collision  terms  of  the  G.B.E. 's,  or  that  those  terms  can  be 
ignored  in  comparison  with  the  collision  terms  of  the  (r-l)-th 
G.B.E.  and  that  N-r  c*  N-r+1. 

To  show  that  Grad's  scheme  satisfies  C2  it  is  sufficient 
to  consider  the  collision  terms.   First  consider  the  second  B.E. 
The  collision  terms  are 

(N-2)y(f^2'3''^123)^23^'^23^^  ^""^   (^-2)/(^l«23'  "  ^123^ ^13^'^13^^ 

Let  us  integrate  with  respect  to  particle  2.   In  the  second  in- 
tegral 2  is  Just  a  parametery and  so  the  order  of  integration 
can  be  reversed.   By  CI,  we  get  (N-2) /(f^i^t  -  f-^^)  V-j^2^^13'^^3  ' 
which,  except  for  the  multiplicative  factor,  is  the  collision 
term  of  the  first  G.B.E.   It  remains  to  show  that  the  other 
integral  gives  zero. 
Note  that 


^i23  =  -2Vaf3-''-V23J^ 


So  we  get 
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/(^2'3'  -^23)%^^23^^^^2  =   "2^l/(^2'^3'  -f2^3)^23^"'23'^^3^^^ 


It  follows  from  the  same  argument  as  is  used  to  find  the 
collisional  invariants  that  each  of  the  three  integrals  on  the 
right  hand  side  of  the  above  equation  vanishes.   The  reason  why 
the  argument  works  is  that  the  integrand  in  each  case  is  symmetric 
in  the  colliding  particles. 

It  is  clear  from  the  above  discussion  that  any  truncation 
which  retains  the  symmetry  in  the  particles  and  satisfies  CI 
satisfies  C2.   This  is  true  of  all  the  truncations  we  have 
discussed,  in  particular  that  of  Grad  in  which 

*i2...r  =  (-1)"(=-1)J   fi+|Z  (-l)^^'[flf2---fjVl...rl^j'  • 

Finally  we  come  to  the  third  condition,  which  we  have 
referred  to  earlier.   The  question  we  are  interested  in  is 
whether  chaos  is  propagated  in  time.   In  other  words,  if  the 
long-range  correlation  is  zero  initially,  will  it  remain  zero 
for  a  macroscopically  long  time?   (If  so,  it  implies  that  the 
r-th  B.E.  together  with  the  initial  condition  gives  the  (r-l)-th 
B.E.  on  integration.)   The  answer  would  be  in  the  affirmative  if 
one  could  prove  an  existence  and  imiqueness  theorem  for  the 
correlation  equation  and  if  the  equation  was  homogeneous  in  the 
correlation.   We  shall  show  that,  under  certain  additional 
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assumptions,  the  condition  of  homogeneitjr  is  satisfied  by  all 

the  three  truncations  we  have  discussed.   Moreover,  the 

correlation  equation  is  the  same  for  the  three  truncations, 

and  it  is  in  this  sense  that  they  are  all  equivalent.   (Strictly 

speaking,  the  correlation  equations  are  not  the  same  in  the 

three  cases,  but  each  of  them  can  be  obtained  from  either  of 

the  other  two  by  using  the  equations  for  lov;er  order  correlations.) 

Recall  that  the  sequence  (7-8)  is  valid  when  e  =  Na-^  <<  1. 
If  the  Mach  number  M  of  the  flow  is  small,  one  can  assume  that 
Ihere  is  a  parameter  5  independent  of  e,  but  depending  on  M 
such  that  the  correlations  are  ordered  in  pov/ers  of  5: 


f(^'^  1 


,(2)    ^2 

(7-12) 


s 

and  so  on.   It  follov;s  from  (7-11)  that  if  the  ^'  are  ordered 

according  to  (7-12),  so  are  the  m-  and  the  A  . 

The  existence  of  statistical  fluctuations  implies  [21] 
that  5  can't  be  smaller  than  £  but  that  5  can  be  >>e.   So,  if 
we  restrict  ourselves  to  those  situations  in  which 


e  <<  6  <<  1  , 
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we  can  use  the  sequence  (7-8)  with  the  ordering  for  the  correla- 
tions.  Also,  for  the  purpose  of  studying  long-range  correlations, 
one  can  ignore  the  force  terms  on   the  left  hand  side  of  the 
equations,  since  it  is  reasonable  to  assume  that  any  two 
specific  molecules  which  are  far  apart  initially  have  a  very 
small  probability  of  colliding  with  each  other  after  a  finite 
amount  of  time.   (The  second  equation  of  (7-8)  with  the  force 
terms  ignored  gives  the  macroscopic  correlation  equations  (6-12), 
but  without  the  potential  parts  of  the  flow.) 

As  an  example,  consider  the  equation  for  ^-.o-s*   ^^  writing 
it  down,  we  shall  put  N-1 ^  N-2  _^  N-3.   Simple  substitution 
and  use  of  the  first  three  equations  of  the  sequence  (7-8) 
gives 


i^j  i^jA 

+  fi(V-#'-^ki'^^^  =  ^V23x+Vl23-'-(^^2,3)} 
+  {a[^^2^3x^^+(l  *^2,3)}-([^-LxV'23]^^+{^123x  +  (l^  2' 3)  3 


Here,  the  brackets  {  ]  mean  that  what  is  enclosed  is  a  collision 
term;  the  underlined  subscripts  refer  to  colliding  particles. 
The  subscript  x  stands  for  the  particle  with  which  each  of  1, 
2,  and  3  collides.   The  symbo]  (1  -^  2,3)  means  similar  terms 
are  to  be  added  vjith  1  replaced  by  2  and  3-   If  nov;  we  ignore 
the  force  terms  and  use  the  ordering,  we  get 
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ignoring  terms  of  order  6  . 

Similarly,  in  the  general  case 

'"^t-   -^fZV    'x!'    -  fV23...rx-^Vl23...r+(l-^2,3,...r) 

r+1 
Ignoring  terms  of  order  6 
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8.    Concluding  Remarks 

We  have  shown  that,  except  in  a  small  part  of  the  configuration 
space,  long-range  macroscopic  correlations  in  a  fluid  are  conserved. 
Since  the  derivation  does  not  depend  on  the  density  of  the  fluid, 
the  equations  of  conservation  are  valid  for  a  liquid  as  well  as 
a  gas.   Zhigulev  and,  independently,  Tsuge  obtained  these 
equations  for  a  rarefied  gas,  and  our  equations  agree  with  theirs 
when  the  terms  containing  the  intermolecular  potential  or  force 
explicitly  are  ignored.   Their  work  as  well  as  ours  shows  that 
contrary  to  the  pessimism  expressed  by  certain  authors  — Monin 
[22],  for  example  —there  is  reason  to  hope  that  a  full-fledged 
kinetic  theory  of  turbulence  can  be  form.ulated.   The  next  step 
tov7ard  that  end  is  to  study  the  properties  of  eqns.  (6-12)-(6-l4) 
and  to  try  to  solve  them.   For  a  rarefied  gas,  Tsuge  has  already 
done  some  of  this  work  by  using  the  two-point  version  of  the  13- 
moment  method.   Finally,  the  results  of  Section  7  show  that  tv/o- 
point,  three-point,  etc.  kinetic  equations  for  a  rarefied  gas 
can  be  obtained  by  generalizing  the  Boltzmajin  equation  in  a 
manner  made  unique  in  a  certain  sense  by  the  imposition  of 
certain  compatibility  criteria. 
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